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We present measurements of the inclusive branching fractions for the decays D^ -^ Xe^Ve and 
D" ^ Xe'^Ue, using 281 pb~* of data collected on the i/)(3770) resonance with the CLEO-c detector. 
We find B{D° -^ Xe+pJ) = (6.46 ± 0.17 ± 0.13)% and 13{D+ -^ Xe+u^) = (16.13 ± 0.20 ± 0.33)%. 
Using the known D meson lifetimes, we obtain the ratio r^+ /^%a ~ 0.985±0.028±0.015, confirming 
isospin invariance at the level of 3%. The positron momentum spectra from D^ and D" have 
consistent shapes. 

PACS numbers: 13.20.Fc,12.38.Qk,14.40.Lb 

The study of inclusive D semileptonic decays is im- tions, one can determine whether there are unobserved 

portant for several reasons. First, by comparing the in- semileptonic decay modes. Previous data suggest that 

elusive branching fractions of the D"*" and D" mesons the lightest vector andpseudoscalar resonances saturate 

with the sum of the measured exclusive branching frac- the hadronic spectra [ill . This may be due to the rela- 



tively low momentum of the daughter s quark, that fa- 
vors the formation of s-wave hadrons. Alternatively, this 
may be an indication that heavy quark effective theory 
may still be valid at the charm quark mass scale [2\. In 
addition, since accurate experimental determinations of 
the D^ and D~^ lifetimes are available [l| , measurements 
of semileptonic branching fractions determine the cor- 
responding semileptonic widths, r^+ and F^o- These 
widths are expected to be equal, modulo small correc- 
tions introduced by electromagnetic effects. Weak anni- 
hilation diagrams can produce more dramatic effects on 
the Cabibbo suppressed partial widths [3| . As these con- 
tributions may also influence the extraction of Vub from 
inclusive B meson semileptonic decays, it is important to 
understand them well. Finally, better knowledge of the 
inclusive positron spectra can be used to improved mod- 
eling of the "cascade" decays 6 — > c ^ se~^i^e and thus is 
important in several measurements of b decays. 

The use of ratios of semileptonic branching fractions 
as a probe of relative lifetimes of the D mesons was sug- 
gested by Pais and Treiman [J]. Indeed the early mea- 
surements of the ratio of the D~^ and 1?° semileptonic 
branching fractions gave the first surprising evidence for 
the lifetime difference between these two charmed mesons 
[a, |6|. Later, the first measurement of the individual 
charged and neutral D inclusive semileptonic branching 
fractions was performed by Mark III [7|, with an overall 
relative error of about 12-16% on the individual branch- 
ing fractions, and 19% on their ratio. The inclusive decay 
D° -^ Xe+Ve was subsequently studied by ARGUS Q 
and CLEO [9| , using the angular correlation between the 
TT"*" emitted in a D*'^ — *■ 7r+I?° decay and the e"*" emitted 
in the subsequent Z?'^ — > Xe^v^ decay. The more precise 
CLEO result has a 5% relative error, dominated by sys- 
tematic uncertainties. Our measurements exploit a clean 
DD sample at threshold, and thus achieves significantly 
smaller systematic errors, through low backgrounds and 
well understood efficiencies. 

We use a 281 pb^^ data sample, collected at the 
i/'(3770) center-of-mass energy {^/s « 3.73 GeV), with 
the CLEO-c detector [lOJ. This detector includes a 
tracking system composed of a six-layer low-mass drift 
chamber and a 47-layer central drift chamber, measur- 
ing charged particle momentum and direction, a state 
of the art CsI(Tl) electromagnetic calorimeter, and a 
Ring Imaging Cherenkov (RICH) hadron identification 
system. All these components are critical to an efficient 
and highly selective electron and positron identification 
algorithm. The charged particle momentum resolution is 
approximately 0.6% at 1 GeV. The CsI(Tl) calorimeter 
measures the electron and photon energies with a reso- 
lution of 2.2% At E = I GeV and 5% at E=im MeV, 
which, combined with the excellent tracking system, pro- 
vides one of the e identification variables, E/p, where 
E is the energy measured in the calorimeter and p is 
the momentum measured in the tracking system. The 



tracking system provides charged particle discrimination 
too, through the measurement of the specific ionization 
dE/dx. Charged particles are also identified over most 
of their momentum range in the RICH detector 11]. In 
particular, RICH identification plays a crucial role at mo- 
menta where the specific ionization bands of two particle 
species cross each other and dE/dx does not provide any 
discrimination power. 

We use a tagging technique similar to the one pioneered 
by the Mark III collaboration [l3|. Details on the tag- 
ging selection procedure are given in Ref . [l3[ . We select 
events containing either the decay D'^ — > K~^tt~ or the 
decay D^ -^ K^n^n". We use only these modes, be- 
cause they have very low background. Note that charge 
conjugate modes are implied throughout this paper. In 
this analysis we exploit the flavor information provided 
by the tagging D: the D~ charge sign provides a flavor 
tag, whereas the charge of the tag daughter K is used for 
D'^ flavor assignment. 

We analyze all the recorded events at the 7/1(3770) 
and retain the events that contain at least one candi- 
date Z?" -^ K~^TT~ or D~ -^ K^TT~TT~ . Two kinematic 
variables are used to select these candidates: the beam- 
constrained mass, Mbc = V-E-boam ~ {^iPiY^ and the en- 
ergy difference Ai?, where /S.E = {T,iEi — -Ebeam), where 
-E'beam represents the beam energy and {Ei,pi) represent 
the 4-vectors of the candidate daughters. For D'^ tags, 
the measured standard deviations {a) in AE and Mbc are 
a{AE) = 6.6 MeV and cr(Mbc) = 1.35 MeV, while for 
D- tags, a{AE) = 5.9 MeV and cr(Mbc) = 1.34 MeV. 
We select events that are within 3 cr of the expected AE 
(0 GeV) and Mbc (M^j) for the channels considered. In 
order to determine the total number of tags, we count 
the events within the selected Ai?-Mbc intervals; then we 
subtract the combinatoric background inferred from two 
3cr sideband regions on both sides of the AE = signal 
peak, with a 2cr gap from the signal interval. The yields 
in the signal region are 48204 D° and 76635 D' . The 
corresponding yields in the sideband region are 788 ± 28 
and 2360 ± 49. We correct the sideband yields with 
scale factors accounting for the relative area of the back- 
ground in the signal and sideband intervals (1.047 for D° 
and 1.23 for D~). The scaling factors are inferred from 
the background component of the Mbc fits. We obtain 
47379 ± 29Z)° tagged events, and 73732 ± 60 Z?" tagged 
events. As we are interested in counting the number of 
signal events, and not in measuring a production rate, 
the errors only reflect the uncertainty in the background 
subtraction. Note that the estimated background is only 
1.7% of the signal for D° and 3.9% for D' . 

For each event selected, we study all the charged tracks 
not used in the tagging mode. We select the ones that are 
well-measured, and whose helical trajectories approach 
the event origin within a distance of 5 mm in the pro- 
jection transverse to the beam and 5 cm in the pro- 
jection along the beam axis. Each track must include 



at least 50% of the hits expected for its momentum. 
Moreover, it must be within the RICH fiducial volume 
(I cos(0)| < 0.8), where 9 is the angle with respect to the 
beams. Finally, we require the charged track momentum 
Ptrack to bc greater than or equal to 0.2 GeV, as the par- 
ticle species separation becomes increasingly difficult at 
low momenta. 

Candidate positrons (and electrons) are selected on the 
basis of a likelihood ratio constructed from three inputs: 
the ratio between the energy deposited in the calorimeter 
and the momentum measured in the tracking system, the 
specific ionization dE/dx measured in the drift chamber, 
and RICH information 1^. Our particle identification 
selection criteria have an average efficiency of 0.95 in the 
momentum region 0.3-1.0 GeV, and 0.71 in the region 
0.2-0.3 GeV. 

The e+ sample contains a small fraction of hadrons 
that pass our selection criteria. As the probability that 
a TT is identified as an e at a given momentum is different 
from the corresponding K to e misidentification proba- 
bility, we need to know the K and tt yields separately to 
subtract this background. We select tt and K samples us- 
ing a particle identification variable (PID) that combines 
RICH a.nd dE/dx information, if the RICH identification 
variable [ll| is available and ptrack > 0.7 GeV; alterna- 
tively PID relies on dE/dx only. The tt sample contains 
also a /i component, as our PID variable is not very se- 
lective; however, as our goal is only to unfold the true e 
spectrum, we do not need to correct for this effect. 

We separate e, tt, and K into "right-sign" and "wrong- 
sign" samples according to their charge correlation to 
the flavor tag. Right-sign assignment is based on the 
expected e charge on the basis of the flavor of the de- 
caying D. The true e populations in the right-sign and 
wrong-sign samples are obtained through an unfolding 
procedure, using the matrix: 
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here n™, n™, n^ represent the raw measured spectra 
in the corresponding particle species, and the coefhcient 
K accounts for the fact that the efficiencies for tt and fi 
selection are not necessarily identical, especially at low 
momenta. The quantities n*, n^, and n^ represent the 
true e, tt and K spectra: the present paper focuses on the 
extraction of n* . As the tt to e misidentification probabil- 
ity is quite small, the effect of a small /i component in the 
measured tt population n™ is negligible. The efficiencies 
Ee, Ctt, and eK account for track finding, track selection 
criteria, and particle identification losses. The tracking 
efficiencies are obtained from a Monte Carlo simulation 
of DD events in the CLEO-c detector. The generator in- 
corporates all the known D decay properties, includes ini- 
tial state radiation (ISR), and final state radiation (FSR) 



effects, the latter are modeled with the program PHO- 
TOS [15|. The particle identification efficiencies are de- 
termined from data. We study the K selection efficiencies 
using a sample of Z?+ -^ K^tt^tt^ decays, and tt selec- 
tion efficiencies using D+ -^ K^tt^tt^ and Kg -^ tt+tt" 
decays. The e+ identification efficiency is extracted from 
a radiative Bhabha sample. A correction for the differ- 
ence between the DD event environment and the simpler 
radiative Bhabha environment (two charged tracks and 
one shower) is derived using a Monte Carlo sample where 
a real electron track stripped from a radiative Bhabha 
event is merged with tracks from a simulated hadronic 
environment. The off-diagonal elements are products of 
tracking efficiencies and particle misidentification proba- 
bilities, where fat is defined as the probability that par- 
ticle b is identified as particle a. The fat parameters are 
determined using e samples from radiative Bhabhas, and 
K and tt from _D+ -^ K^tt^tt^ and K^ -^ tt+tt^. The e 
spectrum from radiative Bhabhas is divided in 50 MeV 
momentum bins to determine the corresponding misiden- 
tification probabilities, whereas the K and tt populations 
are subdivided into 100 MeV momentum bins to reduce 
the statistical uncertainty. The hadron to e^ misidentifi- 
cation probabilities are of the order of 0.1% over most of 
the momentum range and below 1% even in the regions 
where dE/dx separation is less effective. 

There are background sources that are charge symmet- 
ric, mostly produced by tt" Dalitz decays and 7 conver- 
sions. We subtract the wrong sign unfolded yields from 
the corresponding right sign yields to account for them, 
motivated by Monte Carlo studies that confirm the ac- 
curacy of this method. 

In order to subtract the combinatoric background, we 
repeat the unfolding procedure determining the true e"*" 
yields from measured e"*", tt"*", and K'^ samples where 
the tags are selected from /S.E sidebands. The decays 
D^ -^ K^TT' and D^ — > K^tt~tt~ have very little 
background: the 13° sidebands give a combinatoric back- 
ground estimate that is 0.2% of the signal yield and the 
D~ sidebands give a combinatoric background estimate 
that is 1.8% of the signal yield. Table U shows the results 
of the intermediate steps involved in the determination 
of the net e"*" yields. Efficiency corrections increase un- 
folded e+ yields with respect to uncorrected e"*" yields, 
while the subtraction of the contribution from misidenti- 
ficd hadrons reduces them. The former effect is dominant 
for the right-sign positron sample, while it is comparable 
in size to the background subtraction in the wrong-sign 
sample. The final yields, identified as "corrected net e"*"" 
include acceptance corrections related to the (cos < 0.8) 
cut, and doubly-Cabibbo suppressed (DCSD) effects in 
£)" decays. As we are using the charge of the tagging 
D^ , rather than its K daughter charge, this correction 
is not needed in the charged mode. 

In order to extract the partial branching fractions for 
Pe > 0.2 GeV, we evaluate the ratio between the net 



TABLE I: Positron unfolding procedure and corrections. The 
errors reported in the intermediate yields reflect only statis- 
tical uncertainties. 





D+ 


D" 


Signal e+ 






Right-sign 


8275 ± 91 


2239 ±47 


Wrong-sign 


228 ± 15 


233 ± 15 


Right-sign (unfolded) 


9186 ± 103 


2453 ± 54 


Wrong-sign (unfolded) 


231 ± 19 


203 ± 19 


Sideband e+(RS) 


168± 13 


15 ± 4 


Sideband e+(WS) 


11 ±5 


11 ±4 


Net e+ 


8798 ± 105 


2246 ± 57 


Corrected Net e+ 


10998 ± 132 


2827 ± 72 



positron yields corrected for geometric acceptance and 
the net number of tags. While the charge of D~ — > 
K~^TT^7r~ reliably tags the flavor of the charged D, in 
the I?° case the K charge occasionally produces an in- 
correct flavor assignment due to the DCSD D'^ — > K^tt^. 
This effect is estimated on the basis of the known value 
of the parameter tdcsd = NiD^ -^ K+Tr-)/N{D° -^ 
K-TI+) = 0.00362 ±0.00029 J]. 

We have considered several sources of systematic un- 
certainties. There are multiplicative errors that affect 
the overall scale of the spectrum, including tracking ef- 
ficiency or electron identification efficiency, accounting 
for Monte Carlo modeling uncertainties. The uncertain- 
ties in tracking and K and tt identification efficiencies 
are taken from the studies discussed in Ref. [l3(. The 
systematic error on the electron identification efficiency 
(1%) is assessed by comparing radiative Bhabha samples, 
radiative Bhabha tracks embedded in DD Monte Carlo 
samples, and DD Monte Carlo samples. These contri- 
butions are common to D^ and D'^ . In addition, we 
have accounted for the FSR uncertainty by varying its 
amount, with a total systematic error of 0.5%. The last 
multiplicative error is the uncertainty on the number of 
tags, estimated by comparing the number of background 
tags in our signal window from the /S.E sidebands and 
from Mbc sidebands. In addition, there are terms that 
are affected by limited statistics, such as misidentifica- 
tion probabilities, or particle identification efficiencies. 
The systematic uncertainty associated with these terms 
is evaluated with a toy Monte Carlo; we perform 10^ it- 
erations of the unfolding procedure, and vary the matrix 
elements within error. The corresponding relative sys- 
tematic error estimates are 0.56% (statistical errors on 
particle misidentification probability and particle identi- 
fication efficiency) and 0.3% (statistical error on tracking 
efficiency). The uncertainty on the combinatoric back- 
ground, accounted for with the sideband positron sample, 
is negligible compared with these components (< 0.1%), 
because of the excellent purity of the tag samples used. 
Thus the total relative systematic error on the branching 



fraction for pe > 0.2 GeV is 1.7% {D°) and 1.8 % (Z)+). 
The partial branching fractions for pe > 0.2 GeV are 
evaluated as the ratio between the corrected net e yields 
and the net number of tags: 



B{D^ 



Xe" 



B{D" -^ Xe^ 



(14.92 ±0.19,tat± 0.27, 



(5.97 ±0.15,tat± 0.10s 



The yield in the unmeasured region [pe < 0.2 GeV) is 
estimated by fitting the measured spectra with a shape 
derived from Monte Carlo. The semileptonic decays are 
generated with the ISGW form factor model [l6|, with 
parameters tuned to experimental constraints such as 
measured branching fractions, with the procedure de- 
scribed in Ref. 'y\. Final state radiation effects are in- 
cluded in the simulation. We obtain f{pe) = ^^s\{Pe < 
0.2 GeV)/rsi = (7.5 ± 0.5)% for D+ -^ Xe+v^ and 
Arsi(pe < 0.2 GeV)/rsi = (7.7±0.9)% for D° -^ Xe+v^. 
The x^ per degree of freedom is 1.23 for D"*" — > Xe^v^, 
0.75 for 13+ -^ Xe^Ve. We studied the sensitivity of 
our analysis to f{pe) using alternative fitting procedures, 
such as a combination of the dominant exclusive channels 
modeled with different form factors [3|- The fractional 
difference in f{pe) with the various methods considered 
is below 4% and is well within the systematic errors as- 
signed. Note that the relative error in the branching 
fraction introduced by the extrapolation to the unmea- 
sured portion of the spectrum is given S f {pe) / {1 — f (pe)) , 
and thus the systematic error on the total semileptonic 
branching fractions is a about 1%. Upon applying this 
correction, we obtain: 



B{D+ -^ Xe+v^) = (16.13 ± 0.20stat ± 0.33 



sys; 



B{D^ -^ Xe+v,) = (6.46 ± 0.17stat ± 0.13sys)%. 

Using the well-measured lifetimes of the D+ and 
D^ mesons, tjj+ = (1.040 ± 0.007) ps, and T^a = 
(0.4103 ±0.0015) ps [1], we normalize the measured par- 
tial branching fractions to obtain differential semilep- 
tonic widths. Figure [T] shows the differential semileptonic 
widths dT^^/dpe m the laboratory frame, where the D'^ 
momentum is 0.243 GeV, and the D^ momentum is 0.277 
GeV. No final state radiation correction is applied to the 
data points. Table |lT] shows the corresponding numer- 
ical values. The errors shown are evaluated by adding 
the statistical errors and the additive systematic errors 
in quadrature. In addition, an overall multiplicative sys- 
tematic error of about 1.5% needs to be included in de- 
rived quantities such as the total semileptonic width to 
account for overall tracking and particle identification ef- 
ficiency uncertainties. The total inclusive semileptonic 
widths are T{D+ -^ Xe+u^) = 0.1551 ± 0.0020 ± 0.0031 
ps-\ and T{D° -^ Xe+v^) = 0.1574 ± 0.0041 ± 0.0032 
ps^^ . The corresponding ratio of the semileptonic widths 
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FIG. 1: Positron differential semileptonic widths dF^'/dpe for 
the decays D^ — > Xe^Ve, (open squares) and D^ -^ Xe^Ve 
(filled triangles) in the laboratory frame. The errors shown 
include statistical and additive systematic errors. The sym- 
bols for D'^ and D^ spectra are slightly shifted horizontally to 
avoid overlapping. The curves are derived from the fits used 
to extrapolate the measured spectra below the p™'" cut. 



positron differential semileptonic 
"^Mn the laboratory frame. The er- 



TABLE II: D+ and D° 

widths dr/dpe(ps~^GeV~ 

rors shown include statistical errors and additive systematic 

errors . 

^TTGeV) dV/dp^{D+) dV/dp^D") 



0.20 
0.25 
0.30 
0.35 
0.40 
0.45 
0.50 
0.55 
0.60 
0.65 
0.70 
0.75 
0.80 
0.85 
0.90 
0.95 



0.25 
0.30 
0.35 
0.40 
0.45 
0.50 
0.55 
0.60 
0.65 
0.70 
0.75 
0.80 
0.85 
0.90 
0.95 
1.00 



0.1598 ± 
0.2185 ± 
0.2538 ± 
0.2925 ± 
0.3281 ± 
0.3064 ± 
0.3047 ± 
0.2716 ± 
0.2479 ± 
0.1864 ± 
0.1359 ± 
0.0892 ± 
0.0444 ± 
0.0221 ± 
0.0065 ± 
0.0007 ± 



0.0142 
0.0121 
0.0116 
0.0121 
0.0127 
0.0130 
0.0115 
0.0111 
0.0104 
0.0088 
0.0076 
0.0060 
0.0042 
0.0028 
0.0015 
0.0005 



0.1664 
0.1935 
0.2966 
0.3081 
0.3088 
0.3047 
0.2828 
0.2631 
0.2422 
0.1951 
0.1547 
0.0948 
0.0498 
0.0344 
0.0120 
0.0020 



0.0311 
0.0248 
0.0247 
0.0231 
0.0233 
0.0233 
0.0214 
0.0212 
0.0196 
0.0179 
0.0158 
0.0121 
0.0091 
0.0070 
0.0044 
0.0020 



of charged and neutral D mesons is r'^+/r|!)o = 0.985 ± 
0.028 ± 0.015, consistent with isospin invariance. 

Finally, we can compare these widths with the sum of 
the semileptonic decay widths for the pseudoscalar and 
vector hadronic final states recently published by CLEO 
13: BiD+ -^ Xe+z/e)oxci = (15.1 ± 0.5 ± 0.5)% and 



B{D° -> Xe+i/e)exci = (6.1 ± 0.2 ± 0.2)%: the mea- 
sured exclusive modes are consistent with saturating the 
inclusive widths, although there is some room left for 
higher multiplicity modes. The composition of the in- 
clusive hadronic spectra is dominated by the low lying 
resonances in the c —* s and c -^- d, in striking contrast 
with B semileptonic decays, where a sizeable component 
of the inclusive branching fraction is still unaccounted for 

Q. 

In conclusion, we report improved measurements of 
the absolute branching fractions for the inclusive semilep- 



tonic decays B{D~^ 
and B{D° ^ Xe+i 



Xe+i/e) = (16.13 ± 0.20 ± 0.33)% 
= (6.46 ±0.17 ±0.13)%. Using 



the measured D meson lifetimes, the ratio F^+/F^o = 
0.985 ± 0.028 ± 0.015 is extracted, and it is consistent 
with isospin invariance. The shapes of the spectra are 
consistent with one another within error. 



ACKNOWLEDGEMENTS 

We gratefully acknowledge the effort of the CESR staff 
in providing us with excellent luminosity and running 
conditions. This work was supported by the A. P. Sloan 
Foundation, the National Science Foundation, the U.S. 
Department of Energy, and the Natural Sciences and En- 
gineering Research Council of Canada. 



Universitat Bonn, Nussallee 12, D- 



[1] 
[2] 

[s; 

[4] 
[5] 

[6; 
[9; 

[10 

[11 

[12 
[13 
[14 



B 592, 1 (2004). 

Rev. D 52, 2783 (1995) 



Current address 
53115 Bonn 

S. Eidelman et al., Phys. Lett 
D. Scora and N. Isgur, Phys 
arXiv:hep-ph/9503486 . 
S. Bianco, F. L. Fabbri, D. Benson and I. Bigi, Riv. 
Nuovo Cim. 26N7, 1 (2003) arXiv:hep-ex/0309021 . 
A. Pais and S. B. Treiman, Phys. Rev. D 15, 2529 (1977). 
R. H. Schindler et al. [Mark II Collaboration], Phys. Rev. 
D 24, 78 (1981). 

W. Bacino et al. [DELCO Collaboration], Phys. Rev. 
Lett. 45, 329 (1980). 

R. M. Baltrusaitis et al. [Mark III Collaboration], Phys. 
Rev. Lett. 54, 1976 (1985) [Erratum-ibid. 55, 638 
(1985)]. 

H. Albrecht et al. [ARGUS Collaboration], Phys. Lett. B 
374, 249 (1996). 

Y. Kubota et al. [CLEO Collaboration], D 54, 2994 
(1996). 

Y. Kubota et al, Nucl. lustrum. Meth. A 320, 66 (1992). 
M. Art uso et a/., Nucl. lustrum. Meth. A 502, 91 (2003) 
[arXiv:h ep-ex/0209009 . 

J. Adler et al. [Mark III Collaboration], Phys. Rev. Lett. 
62, 1821 (1989). 

Q. He et al. [CLEO Collaboration], Phys. Rev. Lett. 95, 
121801 (2005) arXiv:hep-ex/0504003 . 
T. E. Coan et al. [CLEO Collaboration], Phys. Rev. Lett. 
95, 181802 (2005) ,arXiv:hep-ex/0506052, . 



[15] E. Barberio and Z. Was, Comput. Phys. Commun. 79, [17] G. S. Huang et al. [CLEO Collaboration], Phys. Rev. 

291(1994). Lett. 95, 181801 (2005) ,arXiv:hep-ex/0506053, . 

[16] N. Isgur, D. Scora, B. Grinstein and M. B. Wise, Phys. 

Rev. D 39, 799 (1989). 



